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O ' Abstract 



We study theoretically the thermodynamics, over a broad temperature 
range (5°C to 125 °C), related to hydrated water upon protein unfolding. 



Q , The hydration effect is modeled as interacting dipoles in an external field, 

O ' 

mimicking the influence from the unfolded surfaces on the surrounding wa- 

X 

' ter compared to bulk water. The heat capacity change upon hydration is 

compared with experimental data from Privalov and Makhatadze on four 
different proteins: myoglobin, lysozyme, cytochrome c and ribonuclease. 
Despite the simplicity of the model, it yields good correspondence with 
experiments. With some interest we note that the effective coupling con- 
stants are the same for myoglobin, lysozyme, and cytochrome c, although 
they are slightly different for ribonuclease. 
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1 Introduction 

Proteins consist of 20 different amino acids with a great diversity with regard 
to size, polarity and charge. The understanding of water and interactions with 
water seems to be important in order to understand protein folding in general, 
and the special feature of cold unfolding of several small globular proteins in 
particular [|l]-P]. 

We in the present work will represent the energy difference between the un- 
folded and folded interior, with regard to the water, by mimicking additional hy- 
drogen bonds from which we calculate the hydration heat capacity change upon 
protein unfolding. A justification of the model is the ability for water molecules 
to form an "ice-like" shell ("iceberg" in the terminology of Frank and Evans ||T0| ) 
around apolar surfaces and thus create more hydrogen bonds. Reduction of both 



enthalpy [pT|-[T^ and entropy |T5|, upon apolar hydration seems to be well 



estabhshed [17 . 



However, the protein interior that becomes hydrated upon unfolding, also 
consists of surfaces that has polarity, which means that the surface has permanent 
dipoles and charges. The heat capacity upon purely polar hydration becomes 



2 



surprisingly negative |]T8|,|T^. For apolar surfaces experiments show that the 



hydration contribution to the heat capacity upon solvation is positive. Also for 
proteins where part of the surface is polar this heat capacity is positive. Thus 
for simplicity we in this work will use the apolar "ice-like" shell picture to make 
a effective model for the hydration effect upon protein unfolding. In this way we 
may neglect some crucial features of polar solvation. 

Finally, we apply equilibrium statistical mechanics to the model and calculate 
the hydration heat capacity increment, which we compare with experimental data 



from Privalov and Makhatadze on four different proteins 19 . 



2 Hydration upon protein unfolding 

We will use a refined version of a model first proposed by Hansen et al. |[7|,[20|,pl 



The model studied here was applied by Bakk et al. p2[ on a complete protein 
folding model, but they did not study the hydration effect separately. In this 
work we will study specifically this hydration upon protein unfolding. 

Protein unfolding involves a cavity formation in water with a rearrangement of 



the water molecules surrounding the unfolded protein p3| , [2^ . When estimating 
the solvation energy of exposing the interior of a protein to water, one has to 
calculate the energy difference between hydrated water, associated to the protein, 
and bulk water [|19|. More precisely, the hydration is defined as the transfer of 
a solute from a fixed position in the ideal gas phase to a fixed position in the 



solvent i.e., water in the present case. 

In order to model the effect upon inserting a surface into water, i.e., unfolding 
of a protein, we use the simplified analogy of a classical electrical dipole in an 
external electrical field whose energy is 

Ei = -ecos'd, (1) 

where e is a bending distortion constant. The angle 'd is the polar angle. Eq. (|l|) 



is the hydration model used in the works by Bakk et al. [22,26], and it extends 



the interpretation of the hydration model applied by Hansen et al. [0,^ and 



Bakk et al. 27, 28 



The idea of representing the solvent by dipoles in protein folding was in- 
troduced by Warshel and Levitt [^, and later applications by Russell and 
Warshel ||], Fan et al. and Avbelj [||]. 



In addition to the energy due to the external field [Eq. (|l|)] we will add a 
coupling term 



Ey = 

2 



thus modeling pair interactions between the water molecules, where Jij is the 
coupling constant between water molecules i and j, and Sj is the dipole moment 
of water molecule i. For simplicity we put |sj| = 1. It can be shown that the 
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energy Ei + E2 [Eqs. (|1|) and (0)] in a mean field solution |]33[ can be represented 
as |22|,P11 



E{i3) = Ei + E2 = -(e + bm) cosi? + -bm^, 



(3) 



where bm, with b = Jij, is the mean field coupling between a water molecule 
and its surrounding water molecules, has to be determined self consistently. The 
average dipole moment is defined by m = (cos^). The standard mean field 
solution in an efi^ective field eg = e + bm is 



m = coth (^) ^ 



RT 



(4) 



Note that here and below e and b are energies per mole, as the gas constant R is 
used. 

For dipoles per protein the partition function for the total hydration con- 
tribution upon protein unfolding is 



2'"" f Ei'ff) 

d(f / sirfd exp ' 
Jo 



RT 



n N 



4:nRTsmh{ee/{RT)) [ bm 



exp 



2RT 



N 



(5) 



The total hydration heat capacity change per mole of proteins is 



^C = RT'^^\nZ 



(6) 
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To obtain AC, the self consistent Eq. has to be solved numerically with 
respect to m. Eq. (^) can be obtained from Eq. via 



1 r din Z dee d\nZdm\ 



3 Discussion 

We want to compare the heat capacity change upon unfolding solvation of the 
protein interior with experiments. The proteins considered are myoglobin (Mb), 
lysozyme (Lys), cytochrome c (Cyt), and ribonuclease (Rns) which we compare 
with experimental data from Privalov and Makhatadze on the hydration contri- 
bution to the heat capacity change upon protein unfolding [^. 

The total hydration heat capacity change is shown in Figure |l[ The parameter 
fit to the experimental data agrees quite well with these data. 

Figure |l| shows that the heat capacity has a maximum around 25°C for Mb, 
Lys, and Cyt, while this maximum is shifted to around 50°C for Rns. Also in 
Table I a similar relation is reflected. With some interest we note from Table I 
that both the "electric field" constant e and the coupling constant b are essentially 
the same for all of the four proteins, but there is a small deviation for Rns. 

This small deviation for Rns reflects itself in the ratios between accessible 
polar and total surface area (reported from Makhatadze and Privalov 



AAp/ AAf, which are almost equivalent for Mb, Lys, and Cyt, while this ra- 



Table 1: Parameters, according to Eq. used in Figure |I| for the fitting to the 
experimental hydration data from Privalov and Makhatadze |]19|. The difference 
in water accessible surface area between the unfolded and the folded protein AAj 
is obtained from Makhatadze and Privalov [^. AAp/AAf is the ratio between 
the polar and total accessible surface area. 



protein 


e 


b 


N 


AAt 


AAp/AAt 




(kJmol-i) 


(kJmol-i) 




(A^) 


(%) 


Mb 


2.05 


8.2 


1240 


18250 


36.5 


Lys 


2.05 


8.2 


800 


14090 


39.2 


Cyt 


2.05 


8.2 


740 


11830 


38.2 


Rns 


2.00 


9.0 


500 


13300 


44.8 



tio is significantly larger for Rns. Hence, the parameters e and b may be regarded 
as effective ones for the combined effect of apolar and polar surfaces as discussed 
in Section |l]. 

One notes that Rns differ a bit from the other three proteins considered. This 
may reflect its larger fraction of polar surface which then also can affect quali- 
tative properties. Our model is more like an effective one for a mixed polar and 
apolar surface. Thus features speciflc for polar surfaces are not properly taken 
into account, but are more or less taken into account by adjusting available pa- 
rameters. E.g., our present model may seem to give too small curvature on flgure 
1 for Rns. A reason for this may be the negligence of quantization, which will 
lower the speciflc heat and thus increase its curvature for decreasing tempera- 
tures. The polar (ionic) forces are relatively strong and the hydrogen atom is 
light, which both favor quantum effects. Use of a two-level system |27, 36,37| for 
this kind of problem can thus reflect such quantization. 



4 Conclusion 

We have proposed a thermodynamical model for the hydration of the protein 
interior that becomes exposed to water upon unfolding. To our knowledge this is 
the first model studied and compared to experimental data on the pure hydration 
heat capacity increment over such broad temperature range (5°C to 125 °C). 

Hydration is modeled in an "ice-like" shell analogy, where the water molecules 
are represented by interacting dipoles in an external field. Compared with exper- 



imental data from Privalov and Makhatadze [T^ for the four proteins myoglobin 
(Mb), lysozyme (Lys), cytochrome c (Cyt), and ribonuclease (Rns) the model fits 
quite well. The specific values of the the field coupling constant (e) and dipole 
coupling constant (6) [see Eq. (^)] are the same for Mb, Lys, and Cyt, while it is 
slightly different for Rns. 

In a future model one should take into account qualitative opposite effects, 
with respect to the heat capacity, by apolar and polar hydration. Experimentally 
one finds that the heat capacity change is negative for hydration of purely polar 



surfaces [|1^,[T9|, in contrast to the apolar surfaces where this heat capacity change 
is positive. 

We have reason to expect that the present model, due to its simplicity, also 
may be useful in more complete protein folding models. 
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Figure captions 

Fig. 1. The hydration heat capacity change upon unfolding of four different 



proteins. Experimental data from Privalov and Makhatadze [0. Parameters 



corresponding to Eq. @, for the fit to the experimental data are listed in Table 
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